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ABSTRACT 

We study dynamical evolution of a resonant triple system formed by an inner EMRI 
and an additional outer MBH. The relevant resonant state (A2 — vui ^ const) is supported 
by the relativistic apsidal precession of the inner EMRI, and, unlike standard mean motion 
resonances, the triple system can have a hierarchical orbital configuration (but different from 
the Kozai process). In order to analyze this unusual resonant system, we extend the so-called 
Hamiltonian approach, and derive a mapping from the EMRI-MBH triple system to a simple 
one-dimensional Hamiltonian. With the derived mapping, we make analytical predictions for 
characteristic quantities of the resonance, such as the capture probability, and find that they 
reasonably agree with numerical simulations up to moderate eccentricities. 

Key words: gravitational waves — binaries: close 



1 INTRODUCTION 

In the solar system, orbital resonances are broadly observed at var- 
ious spatial scales (Peale 1986; Murray & Dermott 2000 (here- 
after MD)). For example, Pluto and Neptune have orbital periods 
of 3:2 and their orbital stability is sustained by this simple relation. 
The resonant states with such commensurable orbital periods are 
termed mean motion resonances (MMRs), a nd have been identi- 
fied also among extrasolar planetary sy stems dLissauer et al.|[20lll : 
IPetrovich. Malhotra. & Tremainell2012h . 

In a recent paper (Seto 2012), triple system formed by an 
EMRI and an additional outer massive black hole (MBH) was nu- 
merically studied, using the post-Newtonian (PN) approximation 
(see Fig[T]for the orbital configuration). Here "EMRI" stands for 
"extreme-mass-ratio inspiral" and repre sents an inspiral o f a com- 
pact object (CO) around a MBH (see lOair et aL I ( I2OO4I) for de- 
tail). The numerical simulations were performed mainly from small 
initial eccentricities, and two resonant states were identified with 
A2 — TOi ~ const and 3A2 — vji — 2^1i ~ const. Here A2 repre- 
sents the mean anomaly of the outer MBH around the central MBH. 
The angles zui and f2i are the longitudes of the pericenter and the 
ascending node of the CO. 

Seto (2012) also discussed astronomical aspects for the triples, 
including prospects for gravitational wave and electromagnetic 
wave observations. The expected numbers of resonant captures (not 
the capture probabilities at the resonant encounters) were roughly 
estimated and the mode A2 — ti7i ~ const turned out to occupy the 
majority of the capture events. 

The two resonant states are induced by the relativistic apsi- 
dal precession of the EMRI and do not depend on the inner mean 
anomaly Ai , unlike the standard MMRs in which two terms propor- 



tional to th e inner and outer mean anomalies nearly cancel dPeald 
1 19861 : IMDI) . Consequently, the resonant EMRI-MBH system can 
have a hierarchical orbital configuration and the masses of the two 
MBHs can be comparable. These properties are remarkably differ- 
ent from the standard MMRs where two orbital periods (equiva- 
lently, two semimajor axes) are comparable but the masses of the 
central b ody must be m uch larger than other ones due to orbital 
stability l lGladman|[T993l) . 

In this paper, we focus our analysis to the resonant dynamics 
of the dominant mode A2 — zui ~ const, paying special atten- 



the so-called Hamiltonian a 


pproach in celestial mechanics 1 Sinclaii 


1972 


YodeJl979l:lHenrardl 


19821 iHenrard & Lamaitrell983f Peale 


mt 


MDI). This method has been applied for the standard MMRs. 



Its basic prescription is to extract the essential dynamical degree 
of freedom from the complicated original triple system and map 
the triple system down to a simple one-dimensional system whose 
dynamics is described by a rescaled Hamiltonian (more precisely, 
in a two-dimensional phase space with a canonical variable and 
its conjugate momentum). Our resonance is an unusual mean mo- 
tion resonance, but certainly classified as an eccentricity-type res- 
onance. Therefore, the important dynamical parameters would be 
the inner eccentricity and the resonance angle A2 — rui ~ const. 
Around the resonance, other parameters approximately behave as 
cyclic variables or constants (see e.g. MD). 

So far, various characteristic behaviours of the standard 
MMRs have been successfully explained with the Hamilto- 
nian approach, taking advantage of basic principles on ana- 
lytical mechanics, such as c onservation of adiabatic invariants 
jBorderies & Goldreichlll984|Peale 198 6: MDJ. In this paper, we 
are primarily interested in whether we can suitably extend the 
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Hamiltonian approach for our unusual resonant state A2 — vj-i ~ 
const. If it works well, we can easily make astrophysical arguments 
on the resonant EMRI-MBH systems without using costly numeri- 
cal simulations, and, furthermore, we can better understand the ef- 
ficient analytical approach itself in a perspective different from the 
traditional analyses for the standard MMRs. 

In this paper, by appropriately handling the effects of the rela- 
tivistic apsidal precession, we derive the mapping from the EMRI- 
MBH triple system to the simple Hamiltonians whose forms are 
identical to those used for analyzing the standard MMRs. We then 
make analytical predictions on the dynamical evolution of the hier- 
archical triples around the resonant encounters. We compare these 
predictions with numerical simulations and confirm good agree- 
ments for certain range of the eccentricity e of the inner EMRI. 

This paper is organized as follows. In §2 we summarize basic 
notations, briefly describe our numerical scheme, and provide some 
of representative numerical results around the resonant encounters. 
In §3 we discuss the relativistic apsidal precession. Later, its de- 
pendence on the inner eccentricity e plays a critical role for the 
overall structure of the mapping. In §4, we compare the strengths 
of the first-order term (oc ecos(A2 — ii7i)) and the second-order one 
(oc cos 2(A2 — ^i)) for our resonant state. In §5, we derive the 
mapping mentioned above, by extending the previous studies done 
for the standard MMRs. In the next three sections, using the derived 
mapping, we make analytical predictions on the resonant dynam- 
ics and extensively compare them with numerical simulations. The 
capture rate is examined in §6. In §7, we discuss the gap of the ec- 
centricity observed at a failure of resonant capture. In §8, we study 
resonant encounters for relatively inclined orbits. We summarize 
this paper in §9. 



2 EVOLUTION OF THE SYSTEM 

Our triple system is composed by two MBHs with masses mo, m2 
and a CO of mi(^ mo, m2). The two components mo and mi 
form an inner EMRI and the third one m2 is rotating outside the 
EMRI (see Fig|T]l. For the orbital elements of the triple, we follow 
the positions of mi and m2 relative to the central MBH mo and 
determine the (instantaneous) semimajor axes a; and eccentricities 
ei {I — 1, 2). Since we only handle triples with nearly circular 
outer orbit and the outer eccentricity 62 1) is not important in 
this paper, we put ei — e for simplicity of notation. Except for §8, 
we mainly study coplanar orbital configurations, as shown in Fig[T] 
and define the mean anomalies A; and the longitu des o f pericenters 
zui {I — 1, 2), following the standard convention ( IMDI) . Below, we 
use the geometrical unit G = c — AI = 1 (M = mo + mi + m2: 
the total mass). 

For numerical evolution of the system, we use the three-body 
ADM Hamiltonian Ht b in the post-Newtonian formalism, and ne- 
glect effects of spins. The Hamiltonian is expanded as 



Htb ~ Hn + Hi + H2.5 (1) 



JSchaferl Il987l; Ijaranowski & Schafeil Il997l: 


Lousto & Nakanol 


I2OO8I; Arun et al.l 20091 Galaviz & Braeemann 


2OIOI) (see also 



Moore 1993). Here Hn is the Newtonian term, and Hi is the IPN 
term, namely the leading order relativistic correction. The 2.5PN 
term 7^2.5 is the first dissipative term caused by gravitational radia- 
tion reaction, and invokes the orbital decay of the system. In Eq.llTJ, 
we put the subscript "TB" representing "three-body" to distinguish 
the rescaled Hamiltonian H introduced in §5. 

In the previous paper ( ISetoll2012h , we included the 2PN term 




Figure 1. The coplanar triple system composed by two massive black holes 
(MBHs; and 2) and a compact object (CO; 1). The MBH and the CO 1 
form an inner EMRI. We measure the angular position tui of the pericenter 
of the CO with respect to the fixed reference direction. The outer MBH 
has a nearly circular orbit and its angular position is given by A2 . We can 
similai'ly define the variable Ai for the CO, but it does not appear in our 
resonant vaiiable A2 — roi ^ const. 



case 


mo 


mi 


m2 


f*lin 


8 


d6 
dT 


D 


I 


0.90 





0.1 


50 


0.017 


0.00168 


24.7 


II 


0.98 





0.02 


30 


0.027 


0.0289 


68.0 


III 


0.999 





0.001 


20 


0.035 


1.43 


890 



Table 1. The model parameters adopted in our numerical simulations. The 
initial outer semimajor axis a2in is adjusted to yield a resonant encounter. 
The outer eccentricity is initially set at e2in — and it stays at a small 
value. The left three columns show the basic pai'ameters chai'acterising the 
resonant dynamics, and they are evaluated for coplanar orbits. The ratio 
q/8 is the critical eccentricity for the shift of the resonant angle A2 — roi 
(discussed in §4). The transverse velocity dS/dr and the coefficient D are 
defined in §5, and evaluated with respect to the Hamiltonian i49i for the 
second order resonance. 



H2. But this term is time consuming and less important for our 
resonance. We thus drop it here. 

The equations of motions for the positions xi and momenta 
Pi of the three masses m; {I — 0, 1, 2) are obtained by taking ap- 
propriate partial derivatives of the Hamiltonian. As in Seto (2012), 
we use the new variable s; = Pi/mi to properly handle the motion 
of the CO with mi <^ 1 (including the test particle limit mi — 0). 
These equations are integrated by a Runge-Kutta method with an 
adaptive step size control (Press et al. 1996, and see also Seto & 
Muto20Il for detail). 

In Table 1, we present the model parameters of our numer- 
ical simulations. Since dependence of the resonant dynamics on 
the inner eccentricity e is our central issue, we systematically an- 
alyze it for commonly aixanged sets of parameters such as masses 
(mo,mi,m2) and the initial inner semimajor axis aun. Among 
the three models listed in Table 1, we mainly use models I and II, 
targeting comparable MBHs, and model III is studied for a specific 
purpose in §4. 

In Figs. 2-4, we present samples of typical orbital evolutions 
of model I around the resonant encounters. We set the initial outer 
distance a2in so that the system transverses the resonant condition 
due to the radiational orbital decay. Throughout this paper, we use 
the outer semimajor axis a2 to show the time. This variable 02 is 
monotonically decreasing from its initial value a2in. 

In Fig|2l we show the results from an initial inner eccentricity 
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300 310 320 330 340 350 360 



300 310 320 330 340 350 360 

Figure 2. Evolution of a coplanar EMRI-MBH triple system from an inner 
eccentricity e ~ 0. 1 (model I). The horizontal axis a2 is the semimajor axis 
of the outer MBH m2 and decreases from its initial value a2in = 360 due 
to gravitational wave emission. The three panels show the resonant variable 
A2 — OTi (the upper left one), the inner eccentricity e (the upper right one, in 
modulo 27r) and the inner semimajor axis ai (the bottom one). This system 
encounters the resonance around 02 ~ 340, and the capture is successful. 
The inner eccentricity e starts to increase afterward. 



e,„ ~ 0.1. The test particle mi is resonantly captured by the outer 
MBH binary at the time 02 ~ 340, corresponding to the ratio of 
orbital periods at ~ 17. Incidentally, the inner eccentricity e starts 
to grow and the inner axis ai decreases very slowly. The resonant 
variable A2 — voj soon localizes around ~ 7r/2. 

For the run shown in Fig|3] we set a larger initial eccentricity 
Bin ~ 0.57. The test particle mi is captured into the resonance 
around 02 ~ 280. In contrast to Fig|2] the combination A2 — rui 
now has a large librational amplitude with a small excluded region 
around n. Again, after the resonant capture, the inner eccentricity e 
increases and the semimajor axis decreases. 

In Fig|4l the initial inner eccentricity e^^ ~ 0.57 is close to 
that in FigH] But the initial orbital phases are different between 
Figs. 3 and 4. While evolutions in Figs. 3 and 4 are similar down 
to 02 ~ 280, their subsequent profiles are completely different. 
Around the critical epoch 02 ~ 280, the inner eccentricity e shows 
a large gap in Fig|4l but the EMRI is not captured into the reso- 
nance, as indicated by the rotating variable A2 — ^i- The inner 
semimajor axis ai also has a small gap, but the following Tisserand 
relation (Murray & Dermott 2000, but now for a coplanar system) 
holds nearly smoothly around 02 ~ 280; 



02 
2ai 



+ 



— (1 — e^) ~ const. 
£12 



(2) 



This relation connects the gaps for e and ai in FiglD 

For the standard MMRs, it is well known that the capture 
becomes a stochastic process when we increase the eccentricity 
of the perturbed mass jBor deries & Goldreich 198 4; Peale 1986; 
lMalhotralll988l ; lDermott, Malhotra. & Murravl[l988l;lMDh . In addi- 
tion, the eccentricity shows a gap if the capture is failed. These in- 
teresting characters are successfully explained by the Hamiltonian 
approach. For our unusual resonance, we make detailed analysis on 
these issues later in §6 and 7. 

Our main targets in this paper are the EMRI-MBH triple sys- 
tems in the resonant state A2 — tui ~ const. But it would be 
worth mentioning th at a related resona nt structure was identified 
in the ring of Saturn dPorco etal.lll984 . The Titan ringlet has the 
semimajor axis of 1.29i?s {Rs ~ 6 x lO^km: the radius of Sat- 
urn) and is in the resonant state At — "ojr ~ const with Ti- 





Figure 3. Evolution of a coplanar triple from a large inner eccentricity e ~ 
0.57 (model I). The inner EMRI is resonantly captured by the outer MBH 
around a2 ~ 280. The libration amplitude is large. 
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Figure 4. Evolution of a coplanai' triple from similar orbital parameters 
(e.g. Sin ~ 0.57) as Fig[3]but with a different relative orbital phase (model 
I). The resonant capture is unsuccessful. We can observe gaps of orbital 
parameters at the resonant encounter around a2 ~ 280. 



tan, the largest satellite of Saturn at the distance ~ 20Rs. Here 
ujh is the longitude of the pericenter of the ringlet and Ay is 
the mean anomaly of Titan. The apsidal precession ujr of the 
ringlet is mainly driven by the multiple moments of Saturn (e.g. 
its quadrupole moment; J2 = 1.6 x 10^^). This ringlet has a finite 
eccentricity ej} ~ (2.6 ±1.4) x 10^'* and a radial width ~ 20km. 

Another potentially interesting example is the Earth-Moon 
system. In its early history, the system might be resonantly af- 
fected by the Sun as an outer perturber, through the state Xsun — 
ruMoon ~ co7^st jTouma & WisdonJl998l ; ICuk & Stewar3l2012l) . 



3 RELATIVISTIC APSIDAL PRECESSION 

As demonstrated in the previous section, our resonant state is char- 
acterized by the following relation between the two angular param- 
eters A2 and zui 



77 = A2 — uji ~ const. 

Taking the time derivative of this relation, we have 

712 ~ ri^i 



(3) 



(4) 
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on average (the dot ' representing the time derivative). Here ni is 
the angular frequency of the object mi (I — 1,2) around the central 
MBH mo, and evaluated with Kepler's third law as 



n2 



mo 



1/2 



mo + m2 



1/2 



-1/2 fai\ 



3/2 



(5) 



(6) 



for ?ni ^ mo + m2 = 1. To characterize the hierarchy of the inner 
and outer orbits, we introduce the factor a as 



Q = — < 1. 

a2 



Then the outer frequency is roughly given as 

3/2 

»i2 ~ n\a 



(7) 



(8) 



for mo = 0{1). 

Next we discuss the apsidal precession rate zu\ of the inner 
EMRI. As is well known for Mercury, relativistic correction gener- 
ates the precession with the rate 



3mi 



3/2 



3pni 



(9) 



at the IPN order jLandau & Lifshitzll97ll) . Here, in order to explic 
itly show the relativistic effects, we additionally defined the post' 
Newtonian parameter p of the EMRI as 

_ mo 
ai 



(10) 



In this paper, we only deal with the regime p ^ 1 where the PN 
framework works well. The relativistic precession (|9] depends on 
the eccentricity e as oc (1 — e^)~^. As we see later in §5, this depen- 
dence becomes particularly important for our unusual resonance. 

From Eqs.(|4j(|6} and ([9}, we obtain the following relation for 
the onset of the resonance 



3p 



-1/2 3/2 



1 - e 
or equivalently 

0-2/3 5/3/1 2n2/3 

02 = 3 ' flj^' (1 — e ) ' r 



(11) 



(12) 



The expression for e = was studied in the previous paper 
(Seto 2012, see also Hirata 2011) and we have the relation be- 
tween the PN parameter p and the orbital hierarchy parameter a 



as p ' 



i"^''^/3. For the eccentric cases shown in the previous sec- 
tion, Ea.l lI2b provides 02 — 360 for Fig|2]and a2 — 285 for Figs. 3 
and 4, reasonably reproducing the dependence on the eccentricity 
e. 

Eq.JI lb is obtained by neglecting influence of the distant outer 
MBH m2 and assuming that the precession rate uj\ is dominated 
by the relativistic effect zu\r- Here we evaluate the Newtonian sec- 
ular contribution iJiN due to m2. For moderate eccentricity and 
inclination, the secular effect zuim is estimated as \ML 



o 3/2 

3 m2a, 

^l'^ = 4 3 1/2 - 



Then, at the critical distance ( I12t . we have 

zu\N 9 mom2 
^ 4" 



< 1 



(13) 



(14) 



VJlr 4 ai 

with p = mo/ffli <C 1 and m2 < 1. Therefore, the Newtonian 
contribution for the precession zu-i would be much smaller than the 



relativistic one. The distant outer body m2 also has a IPN effect 
for the precession vji (see the IPN interaction term in Naoz et al. 
2012). But its magnitude is 0(a^'^) times smaller than Eq.([9]l, and 
not important for the precession zu\. We hereafter put 



as already assumed to derive Eg. ( 112) . 



(15) 



4 COMPARISON BETWEEN THE FIRST AND SECOND 
ORDER RESONANCES 

The gravitational interaction between the inner and outer orbits of 
a triple system has been perturbatively analyze with the disturbing 
function ( IMDI) . For our resonant state 77 ~ const in a coplanar 
configuration, the relevant element of the disturbing function is ex- 
panded as 

C = e^Ci,o cos ?7 + 6^6*2,0 cos 2j7, (16) 

where we take the terms up to the order O(e^). The functions Ci,o 
and C2,() depend on the hierarchy parameter a ^ 1 of the orbital 
configuration. They are explicitly given as 



2 / ada , i\ ,(0) . ^ 15 3 



(17) 



„ , 3 , Zada , 3q da 
= I 4 + ^ + 



^1/2 



15 



105 



8 



64 



with the Laplace coefficients h^p {a). In Eg.llITt the first term 2a/3 
is the indirect part and is canceled by the 0(a) term of its direct 
part. As a result, the function Ci,o oc o? has a stronger dependence 
on the parameter ot{<^ 1) than the counterpart 6*2,0 c< ot^ . Actu- 
ally, the second-order one C2,o has the lowest power of a among 
the resonant terms in the form cos Nr) with N ^ 1. We hereafter 
neglect the term oc in Eq.JISt and put 



. 15 3 

C — OL e 

^ 16 



cos 77 — 2— cos 2r) 

a 



(19) 



This expression shows that the second order term can dominate 
the first order one even at a small eccentricity e ~ a, due to the 
hierarchy of the system a ^ 1. 

Interestingly, the competition of the two terms can be directly 
observed as a shift of the mean angle rj of libration, during the 
resonant amplification of the inner eccentricity e. We now discuss 
this in some detail. For simplicity, we assume that the dissipative 
evolution is negligible during one libration period. 

First, the system around the resonance can be effectively re- 
duced to one dimensional system (with the variable oc 77 and its 
conjugate momentum oc e^, see §5 for detail). The effective Hamil- 
tonian has the resonant term oc and the variable 77 appears only 
in this term. Then, from the canonical equation, we should have 



— oc ^ = 
dt drj 



(20) 



at the equilibrium point {ee,Tje). Thus, for given equilibrium value 
e = ee, we associate the corresponding equilibrium angle 77^ as the 
minimum of the following potential V{(x Q 



2e 

V = — cos 77 H cos 277. 

a 



(21) 



The shape of this potential is shown in Fig|5]for representative val- 
ues of the ratio ee/a. The positions of the potential minima qualita- 
tively change at the critical value Ce = q/8. In Table 1, we present 
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its value for models I-III. At < a/8, the potential V is domi- 
nated by the first order term and we have the equilibrium angle 



»7e = 0, (ee < a/8). 



(22) 



When increasing beyond the critical value a/8, the angle ije 
starts to move as 



rje = ± arccos(a/8ee), (ce > a/8). 



(23) 



We have rje — ±7r/2 for Ce 3> a/8, dominated by the second 
order term in Ea.( l2U . 

Now we examine our simple model i22\ and i23\ for the equi- 
librium resonant angle, by using numerical simulations. In Fig|6l 
we show the evolution of orbital parameters for model III. Owing 
to its small outer mass m2, the forced eccentricity is small at the 
early stage, and this model allows us to make a suitable demonstra- 
tion for the present analysis. 

The EMRI is resonantly captured by the outer MBH binary 
around a2 ~ 55, and its eccentricity e starts to increase afterward. 
Here the critical eccentricity for the onset of the shift of the equilib- 
rium angle is a/8 ~ 0.035. Since the libration can be regarded as 
a circulation around the equilibrium point (ce, rje), the mean value 
of the libration would be close to the equilibrium point (ee, ry^), at 
least for a small libration amplitude. To directly show the shift of 
the angle r^e thorough the resonant amplification of the inner ec- 
centricity e, we plot the combination (e, Tj) in Fig |7] together with 
the analytical model {12\ and ( |23| l. 

Fig[7]shows that, even though the dissipative time scale is not 
sufficiently long compared with the libration period, the simple an- 
alytical prediction shows a good agreement with the numerical one. 
For a larger libration amplitude as in Fig|3] the potential wall of V 
around 77 = (see Fig|5]( is easily crossed over, and the angle 77 
moves around a broad region, leaving a small excluded regime near 

?7 = TT. 

Note that also in Fig[2l the angle 77 finally localizes around 
7r/2, as in the case of Fig.6. But this should be regarded as a mere 
coincidence. Later in §6 and 7, we deal with a large sample of nu- 
merical simulations for models I and II. Among them, there are no 
definite asymmetries for the preference of the two potential minima 
at 77 = 7r/2 and — 7r/2 (equivalently 37r/2). 

In summary, due to the hierarchy of the orbits with a <c; 1, the 
second order term oc e^a^ cos 277 could become more important 
than the first order one oc ea'^ cos 77, even for a small eccentricity 
e ~ a. We can observe the resultant shift of the mean (equilibrium) 
angle 77e during the resonant amplification of the inner eccentricity 
e. 



5 HAMILTONIAN APPROACH 

In this section, we apply the Hamiltonian approach for the reso- 
nant dynamics of the EMRI-MBH triple systems with A2 — roi ~ 
const. By taking appropriate set of conjugate variables, the dy- 
namics around the resonant encounter can be reduced t o a sim- 
ple one dimensional system (ISinclaiilll972l [Yoder 1 19791 : iHenrardI 



1 1 9821 : iHenrard & LamaitrdfigSSl: |Pealelll986t ImdI) . For the stan- 
dard MMRs such as 2:1 or 3:1 resonances, this approach success- 
fully exp lains characteristic phenomenon around the res onant en- 
counters jBorderies & Goldreichiri984|Pealelll986l:lMDl) . Our aim 
here is to extend it for our unusual resonance. As detailed descrip- 
tions of the approach for the standard MMRs can be found in the 
literature and many of them are shared with our resonance, it would 
be unfruitful to lengthily expound all the involved steps. We rather 




Figure 5. The effective potential V defined in Eq.(2T). Its shapes are plot- 
ted for representative values of Ce/o. The equilibrium angle 77^ would be 
the minimum of the potential For Ce/o < 1/8, the potential has the sin- 
gle minimum at r; = 0, dominated by the first order terni. In contrast, for 
Ge/a > 1/8, we have two minima at 77 ^ 0, reflecting the second order 
component. For larger ee/a, the minimal points approach ±77/2. 
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h 
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Figure 6. Evolution of the orbital elements of the coplanar inner EMRI 
for model III. Due to the small outer mass 7712 , the forced components can 
be suppressed, compai'ed with cases I and II. We have a small initial inner 
eccentricity e ^ 0.0065. The EMRI is captured by the outer MBH binaiy 
around 02 ~ 55. We can clearly observe the shift of the equilibrium angle 



C 

■3 

as 



II 




0.20 



Figure 7. CoiTespondence between the inner eccentricity e and the reso- 
nant angle = A2 — roi. The points are obtained from the ran shown 
in Fig.6. The red curves are the analytical predictions \22) and (23) for the 
equilibrium points (ee , r)e) with the transition eccentricity at q/8 = 0.035. 
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follow the comprehensive formulation given in §8.8 of lMDl and ex- 
plain the modifications necessary for our specific resonance. 



5.1 Simplified Hamiltonian for Second Order Resonance 

Based on the results in the previous section, we analyze the second 
order resonance with the resonant variable 



(24) 



identical to Eq.(8.78) of MD with k — 2. While we are mainly 
interested in the specific case j = 2, we do not fix the parameter 
j at this stage, in order to enable a simple comparison with the 
standard second-order MMRs corresponding to j > 2. 

As explained in MD, the variable 6i has the conjugate mo- 
mentum Qi defined by 



mi 
mi 



Vmooi I 1 



Vrnooi 



(25) 
(26) 



We should notice that this momentum is directly related to the inner 
eccentricity e as Si oc e^. In our analytical studies below, we make 
perturbative expansions, assuming e ^ 1. 

Among multiple terms in the Hamiltonian (8.98) of MD (de- 
noted as H), the key element for our unusual resonance is the fol- 
lowing one 



kQiujs 



X 



(27) 



with k = 2 for the present analysis. Here, the notation vjsec in 
MD represents the secular precession rate of the inner peric enter 
and is identical to the relativistic apsidal precession rbir under our 
prescription in §3 (hereafter using zuir in stead of uisec)- 
With respect to the canonical equation 



d6>i 
'dt 



m 

dOi' 



(28) 



the term X in the total Hamiltonian H has a role to provide the sec- 
ular contribution —kvoir for the time derivative d^i /dt. Therefore, 
we should have the equation below 



dX 



-kiijir 



(29) 



Meanwhile, as given in Eq.([9l(, the relativistic precession rate rbir 
has the following form at IPN order 



UJir — 



3pni(l + e ), 



(30) 



and the rate vjir itself depends on the conjugate momentum Qi oc 
e^. Thus we have the following perturbative solution X for EQ.( l29t 



X 



'3kpniQi 1 



-SkpriiQi 



261 



(31) 



(32) 



expanded in terms of the momentum B, instead of the eccentricity 
e. Note that this solution is different from the naive expression l |27) 
that is perturbatively expanded as 



kQiujsec = — 3fcpni6i 1 + 



461 



(33) 



The quadratic term oc Of plays a critical role for our resonance, as 
we see in the next subsection. This term originates from the depen- 
dence VDir OC (1 — e^)^^. 

One might has an impression that the present derivation for 
Eg. (132b is phenomenological, as it is constructed to reproduce the 
desired precession rate rbir- But we can actually derive the term 
(proportional to + 3e*/4 + ■ • ■) identical to X in Eq.ll32t. start- 
ing directly from the IPN Hamiltonian Hi in Eq.Q (see Eq.(28) 
in Naoz et al. 2012). We took the above route to elucidate the mod- 
ification relative to the typical analysis for the standard MMRs. 

With the explicit form of the relativistic correction X in hand, 
we can next apply the standard arguments in MD to derive a sim- 
plified Hamiltonian for MMRs. After some calculations (e.g. intro- 
ducing the new conjugate variables O'l = 6'i/2 and F = 20i), 
we have the following Hamiltonian (corresponding to Eq.(8.102) 
of MD) 



ar + ;3r^ + 2ercos 29'i. 



Here the coefficients a, jS and e are given as 
_ _ ( j — 2)ni — jn2 + 6pni 



C,. 



(i-2f 
mia\ 



+ 



J 



+ 3 



2 ' 

miflj 



m2 

-2^^l Q. 

mo 



(34) 



(35) 



(36) 



(37) 



In Eqs.ll35ll and ( 136b . the terms proportional to the FN parameter p 
clearly show the relativistic corrections. The factor = C2.0 

for j = 2 was already given in Eg. (118b. In the right-hand side of 
Eg. (136b. the first parenthesis appears in the standard MMRs and has 
its origin in the Keplarian terms in the triple system (see MD). Its 
second term (oc p) is due to the quadratic term oc 0f in the secular 
correction X for the relativistic apsidal precession. 

We further make transformation of variables as follows 



El 

2e 



ei/3 



r = 2it 



(38) 



(39) 



(40) 



(41) 



(42) 



(O'l (e<0) 
^ \9'i+n (e > 0) ' 

and finally obtain the rescaled Hamiltonian 

= <I>2 + 5$ + $cos2(^ 
with the single parameter 5 defined by 

21 

The associated canonical equations are written as 

rf$ _ dH d(j> _ dH 
dr d(j) ' dr 9$ ' 

The rescaled Hamiltonian ( 140b is slightly different f rom the 
related exp ression (8.116) in MP, but identical to those in lOuillerj 
( l2006h and lMustill & Wvatll ( l20Ilh . We adopt the present form, in 
order to use these two references later and discuss whether evolu- 
tion of the parameter 5 can be regarded as adiabatic for our resonant 
dynamics. 

Roughly speaking, this parameter S represents an effective dis- 
tance to the resonance. Due to the GW emission, the orbits of the 
EMRI-MBH triple system decay gradually, and the parameter S 
varies accordingly. 

We now estimate the transition rate dS / dr. First, apart shortly 
from the triple systems, we consider a simple binary with a semi- 
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major axis a, an eccentricity e and mass es m, m! . Its orbital decay 
rate da/dt by GW emission is given as ( |Petersll964 



da 64: mm' {m + m') 
It ^ 5" a3(l - e2)7/2 



, 73 2 37 4 
^ + 24^ + 96^ 



(43) 



Next, for our triples, we assume that, before the resonant en- 
counters, the EMRI and MBH binary independently evolve with 
Eq.(|43b. Then we obtain 



da 



/2 -13/2 {(i - 2)ai + lOmo} 



mimQ flj 

^ ^ 73 2 ^ 37 4 
^+24^ +96^ 



(l_e2)7/2 



' J JTlo 771202 



11/2 



(44) 



for 62 = 0. With the scaled time r, we can obtain the transition rate 
as 



(45) 



dS 1 da dt 1 da 
dr ^ 2^'dtdT ^ 4^'dt' 

For a given EMRI-MBH triple around the resonant encounter, 
we can now analyze its evolution through the one-dimensional 
rescaled Hamiltonian J40t . The information of the original triple 
system is converted to (i) the new variables (0, (ii) the param- 
eter S and (iii) its time derivative dS/dr. In practice, this mapping 
can be made with Eqs.l|35ll-l|39t and ll44t-l|45t. In the next sub- 
section, we concretely study the relation in the test particle limit 
mi — >■ 0. But, here, we derive a result valid also for mi 7^ 0. 

To realize a capture (i.e. transition of from rotation to libra- 
tion) with Eq.ll42t. the resonance should be crossed in the direction 
dS/dr < (Peal 1986; MD). In the cases of standard MMRs, 
this corresponds to relatively approaching orbits. For example, to 
be captured into the 3:2 resonance, the ratio of the orbital periods 
should change in the direction of 1.6 — 1.5 not 1.4 — >■ 1.5. With 
Eqs.JI2t and l l44b for j = 2, the inequality dS/dr < is rewritten 



ai 
a2 



5/2(i + Z^e2 + 
24 



37 
96^ 



/ mpmi \ 
V m2 / 



1/2 



.(46) 



For the specific case e = 0, this expression agrees with that derived 
and examined in Seto (2012). Note that our labels (0, 1, 2) for the 
three masses are different from those in Seto (2012). 



5.2 Test Particle Limit 

Here we discuss the mapping between the EMRI-MBH triple sys- 
tem and the simplified Hamiltonian system j40| ) in the test particle 
limit mi 0. In this limit, we can easily control the relative or- 
bital evolutions of the triple system in numerical simulations, and, 
furthermore, the role of the post-Newtonian corrections becomes 
transparent. 

From Eqs.( l26t ( l36t ( l37t and ( I38t . the inner eccentricity e is re- 
lated to the momentum <I> as 



2 

mofli 
20m2a'' 



U - 2)^ 



m2a2 



8p 



(47) 



In this relation, we pay our attention to the dependence of the mass 
parameter mi. We can put p = in the traditional analysis of the 
standard MMRs with j > 2 (see Eq.(8.109) of MD). However, for 
our unusual one with j = 2, the mapping l l47t becomes singular 
$ oc mie^ in the limit mi — >■ 0, if the relativistic effect is dropped 
with p = 0. Therefore, interestingly, the regularity of the mapping 
( 147 1 is maintained by the post-Newtonian correction (oc p) for our 
resonance with j = 2 as 



2pmoa2e'^ 
5m2af 



(48) 



without depending on mi. As mentioned earlier, the post- 
Newtonian term in Eq.(l47t comes from the quadratic term of the 
momentum $ in Eq.i32t and intrinsically from the dependence of 
the precession rate on the eccentricity as shown in Eq.(|9). 

Now we derive formulae specifically for j — 2 with the 
rescaled Hamiltonian 



H = +(5$ + $cos20. 



(49) 



The variable is related to the original resonant angle 77 = A2 — taJ2 
as 

= 7) + TT. (50) 

After some algebra with Eq.( II2) . the principal quantities for the 
rescaled Hamiltonian are given by the original parameters as 

2 



D ; 



d5_ 

dr 



2ai 



45mom2 



512m: 



3/2 



3i/3375m2a; 



13/6 ■ 



(51) 
(52) 

(53) 



Note that the transit speed dS/dr is nearly constant around the res- 
onant encounter, and we omit the expression for the time dependent 
parameter 5 itself. 

In Table. I, we provide the transit speed dS/dr as well as the 
coefficient D. The former is an useful measure to discuss the adi- 
abaticity of the time evolution of the parameter 5 at the resonant 
encounter 

For comparison, including only the first order resonant term 
oc e in Eq.(|19b, we derive the relevant expressions for the test par- 
ticle limit. In Appendix A, we summarize the results. Again, we 
have a regular mapping between the momentum $ and the inner 
eccentricity e, due to the PN correction. 

For a Newtonian apsidal precession induced by multiple mo- 
ments of masses, the precession rate at e < C 1 generally has correc- 
tion for the eccentricity e starting from ( IStemell939f) . Therefore, 
the mapping to a corresponding rescaled Hamiltonian becomes reg- 
ular, as for the relativistic one discussed above. For example, the 
precession rate of a test particle due to the quadrupole moment J2 
of the central body is expanded as ri? oc J2(l + 2e^ + •••). 

In the next three sections, using the mapping from the EMRI- 
MBH triple system, we make quantitative predictions on the reso- 
nant dynamics and compare them with numerical simulations. Be- 
low, we limit our analysis to the test particle limit mi — s> 0. 



6 CAPTURE RATE 

In this section, we study whether the analytical model based on the 
rescaled Hamiltonian can reproduce the capture rate estimated from 
numerical simulations. 

For this comparison, we obtained the capture rate from the 
numerical side in the following manner First, for models I and II, 
we took various (~ 13) inner initial eccentricities ei„ between ~ 
and ~ 0.7. For each eccentricity, we assigned an initial outer radius 
a2in larger than Eq.( ll2) to assure a resonant encounter, and made 
20 runs starting from randomly distributed relative orbital phases. 
Therefore, the total number of the runs is ~ 2 x 13 x 20 ~ 500. We 
judge a run as a resonantly captured event when the angle rj = X2 — 
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Figure 8. Resonant capture rate of an inner EMRI by a coplanar outer 
MBH. The results are given as functions of tlie inner eccentricity for 
models I and II. The solid curves are the analytical predictions P2 for the 
second order term oc cos 2r] and the dashed curves are Pi for the first- 
order term cx: cos The circles are numerical results for minimum dura- 
tion 02/10, and the additional open squares are for 02/20 (shown only if 
different). Each of them is estimated from 20 runs. 

Tui has a single excluded range A77 larger than tt/IO for minimum 
duration 02/10 (and also 02/20 for comparison) in terms of the 
decaying outer radius 02 (see e.g. Figs. 2 and 3). An MMR is often 
identified with a sharp concentration of a resonance angle, as found 
in Fig|2] But, here, we employed the criteria for At; to handle the 
resonant state with large libration amplitudes, as demonstrated in 
FiglD 

Then, for each initial eccentricity ei„, we counted the number 
of captured events among the associated 20 runs and roughly ob- 
tained the capture rate. We provide the numerical results in Fig. 8 
with the blue circles for minimum duration 02/10. Up to moderate 
eccentricity, we obtain the same results for the minimum duration 
02/20. Only when they are different, we added the latter with the 
open squares. . As shown in Figs. 2-4, the eccentricities of the EM- 
RIs are always oscillating to some extent. To handle this, we took 
time averaged eccentricity for each run at its early stage, and sub- 
sequently evaluated the mean value among the 20 runs. The initial 
eccentricities in Fig. 8 are made up in this way. 

Next we analytically estimate the capture rate through the 
rescaled Hamiltonian l |49t , following Borderies & Goldreich 
(1984). In Appendix B, we briefly describe their results. As we 
have already discussed, the second order resonance is relevant for 
our systems (unless Bin is less than 0(10^^), see Table 1). There- 
fore, we mainly use their results for k = 2 (see §B1). 

For an initial eccentricity ei„, the analytical rate P2 is obtained 
through the projected initial momentum 

$„, = DeL (54) 
for the rescaled Hamiltonian. There is a critical value <!>„■ — 1, 



7 

Figure 9. Gaps of the eccentricities e of inner EMRIs with failed captures 
(as shown in Fig|4). With the solid curves, we plot the analytic correspon- 
dences between the gapped eccentricities before (e^^) and after (eji„) the 
resonant encounter for models I and II. The circles are the mean eccentrici- 
ties of the numerical runs accompanied by the gaps. 



comsponding to ei„ — ^ ^cr/ D — 0.20 (model 1) and 0.12 
(model II) that are much larger than a/8 — 0(10^^) (see Table 
1). The capture occurs at 100% for < ^cr- Meanwhile, for 
$171 > $c?-, the capture becomes a stochastic process with the rate 
P2 defined implicitly by Eas.( IB2t and ( |B3l l. 

In Fig. 8, with the solid curves, we present the analytical rates 
P2. At Ein < 0.4, they show reasonable agreements with the nu- 
merical ones. But, at larger eccentricities ein > 0.5, we have sig- 
nificant discrepancies. This is not surprising, since we made, at 
least, various approximations, valid only for <^ 1. For refer- 
ence, we also show the rate Pi expected for the fist order resonance 
(see §B.2), but it poorly fits the numerical data, as expected. Note 
also that, for higher eccentricities, the numerical results are affected 
by the applied conditions for identification of the resonances. 

As we explained earlier, the critical eccentricity ^cr/D 
characterises the overall shape of the capture rate. Here we should 
comment on its magnitude for the standard second-order MMRs 
(with the variable J24t at j > 2). For these resonances, dynami- 
cal stability of orbits requires m2 <Si mo (assuming mi < m2). 
Then we can show a scaling behaviour D = $/e^ oc l/m2 and 
obtain the critical eccentricity $cr/D <C 1 much smaller than 
our hierarchical one with mo/m2 = 0{1) (as in model I). There- 
fore, for the standard second-order MMRs, perturbative expansion 
of the eccentricity is more effective in the regime where the reso- 
nant capture is probable (e.g. P2 > 0.1), unlike our hierarchical 
one with larger -y/ $cr/D. We can make similar arguments for the 
first-order resonances. 

For the analytical predictions in Appendix B, we fully use the 
arguments based on the adiabatic invariant that is conserved for 
a transit speed \dS/dT\ m uch smaller than the libration frequency 
jLandau & Lifshitzll 19691) . To examine the impacts of finiteness of 
\dS/dr\ on the resonant dynamics, Quillen (2006) and Mustill & 
Wyatt (201 1) numerically studied dependence of the capture rate on 
the transition speed | di/dr | . Their results (see e. g. Fig. 2 in Mustill 
& Wyatt 201 1) indicate that the adiabatic approximation would be 
efficient for |d(5/dr| < 0.1. As shown in Table. 1, two models I 
and II well satisfy this criteria. For a coplanar EMRI-MBH triple 
of comparable MBHs (m2/mo = 0(1)) with converging orbits 
dS/dr < 0, we generally have \d5/dT\ <^ 0.1, unless the target 
EMRI is highly relativistic. 
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7 GAP OF INNER ECCENTRICITY 

It is well known that, for the standard MMRs, the eccentricity 
of a perturbed body shows a gap when the resonance is encoun- 
tered but capture results in failure (e.g. Peale 1986; MD; see also 
Amaro-Seoane et al. 2012). This phenomena is well explained by 
the Hamiltonian approach in the associated phase space, as a rapid 
change of rotational motion at the separatrix crossing. In Fig|4] we 
can observe a similar gap of the eccentricity for the run without a 
resonant capture. 

In order to further examine validity of our Hamiltonian model 
extended for the unusual resonance, we analyze the correspondence 
of the two eccentricities; ei„ (before the encounter) and e/i„ (af- 
ter the encounter). We derive an analytical prediction using the 
rescaled second-order Hamiltonian ( 149) and compare it with the 
numerical simulations done in §6. 

We first discuss the analytical approach in which the corre- 
spondence between the two eccentricities is equivalent to the rela- 
tion between the initial momentum $i„ ~ Del„ and the final one 
"J* fin = ^C/in both far from the resonant encounter. For a given 
initial momentum the parameter 5 at the separatrix crossing 
is given in a somewhat complicated form as in Eqs. dBSt and iB4t . 
But, because of a simple expression for a define integral , we have 
the fo llowing concise relation between the two momenta jMalhotral 
Il988h 



fin 



(55) 



Note that the separatrix relevant for our analysis is formed at 5 < 

— 1 where the capture rate becomes less than unity. In Eg. (155), we 
have (^in ~ 1 and $ fin = for 5 = —1. 

We can now obtain the desired correspondence — s> "I> fin 
(equivalently ei„ — ^ e/i„) through the intervening parameter S{< 

— 1). In Fig.9, we show this analytical correspondence for models 
I and II with the solid curves. 

We also analyze the samples of the numerical simulations de- 
scribed in the previous section. In Fig.9, the numerical results are 
presented with the filled circles. In numerical data, formation of 
a gap can be easily identified as a instantaneous event, compared 
with continuation of a resonant state. We can observe small sys- 
tematic deviations between the analytical and numerical results. 
But, as a whole, the simple analytical predictions show reasonable 
agreements with the numerical ones that were obtained after rather 
complicated dynamical evolutions. 



8 INCLINED ORBITS 

So far we have discussed the resonant dynamics for coplanar orbits. 
In this section, we extend our study to inclined orbits. We use the 
parameter i as the relative inclination angle between the inner and 
outer orbits. 

In Fig. 10 we provide a numerical example for inclined orbital 
configurations. This triple system has a small initial inclination an- 
gle i — 5.8°, but its initial semimajor axes ai, a2 and eccentricity 
e are close to those in Fig.2. We find that the overall evolution of 
the three quantities A2 — roi, e and ai are similar to Fig.2. 

Note that the inclination angle i stays nearly at a constant 
value. This is in accord with the simplified Hamiltonian approach 
that has only two dynamically important variables e and A2 — taJi 
in the present eccentricity resonance. 

In Fig. 1 1, we show the results from a larger inclination angle 
i = 60°. In the lower right panel, evolution of the inclination angle 
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Figure 10. Evolution of an EMRI-MBH triple system from an inner eccen- 
tricity e = 0.10 (model I). This system has initial parameters similar to 
Fig.2, but with an initial relative inclination i ~ 5.8°. The EMRI is cap- 
tured into resonance around 02 ~ 340. The inclination angle i is within 
5.7° < i < 6.0°. 
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Figure 11. Evolution of an EMRI-MBH triple system from an inner eccen- 
tricity e 0.39 and relative inclination i ^ 60° (model I). This system 
encounters the resonance around 02 ~ 305 where the inner eccentricity e 
shows larger variation Ae ~ 0.2 compared with A cos i ~ 0.03. 



i is presented in a geometric form cos i. We can observe oscillation 
of cos i. But, around the resonant encounter 02 ~ 305, its ampli- 
tude is much smaller than that of the eccentricity e. 

In the analytical Hamiltonian approach, we need to clarify 
how the resonant interaction depends on the inclination angle. 
Here, based on the above numerical demonstrations, we make an 
approximation that the inclination angle i is constant around the 
resonant encounter. In the previous case for coplanar orbits, the 
disturbing function has the following second-order resonance term 



15a^e^ 
C = — z — ^'^^ ^""l- 



(56) 



Here we neglected subleading contributions of o(a^e^). In celestial 
mechanics, the effects of the inclination on the disturbing function 
are often handled perturbatively with the expansion parameter s = 
sin[i/2]. But, here, we are interested in highly inclined orbits with 
s = 0(1), well beyond the perturbative regime s ^ 1. 

We should notice that the factor e for the resonant term in 
the Hamiltonian J34t is contributed by all the terms proportional to 
cos 2ri among the disturbing function. Fortunately, we can readily 
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collect the terms at the lowest order a e cos 277 as follows; 



-(1 



2s^ + s^) cos 2?7 



8 



s*[i/2] cos 2r;. 



(57) 
(58) 



Therefore, with respect to the original Hamiltonian given 
in Ea.J34t and at the order 0{o?e^) of the resonant interaction, we 
just need to multiply the factor cos^(i/2) to the parameter e that 
was defined in Ea.J37t for the coplanar system. Under the present 
approximation i = const, this is basically what we need to do 
for dealing with the inclined orbits. Accordingly, for the rescaled 
Hamiltonian {49} in the test particle limit (see §5.2), the coefficients 
D in Eq.l[52) and the transition speed dS/dr are given as 

2ai 



D = 
_ 

dr ~ 



45mom2 cos'*[i/2] ' 



512m: 



3/2 



3i/3375m,2af/''cos8[i/2] 



(59) 



(60) 



With these expressions, it is straightforward to apply the previous 
analytical methods in §6 and 7 to inclined orbits. 

Now, we compare these analytical predictions with numerical 
simulations for triple systems. Below, we fix the inclination angle 
at a relatively large value i — 60° and perfomed a large number 
of simulations for two models I and II. Even with the strong de- 
pendence on the inclination oc l/cos**[i/2] = (4/3)* — 3.2, the 
transition speeds dS/dr are less than O.I both with models I and II 
(see Table 1), and the adiabatic approximation would be still effec- 
tive for analyzing resonant dynamics. 

Note also that due to the relativistic apsidal precession, the 
Kozai process ('Kozai"l962'; 'Lido \i Il962h does not work here 
(jHolman, Touma, & Tremaine 1997: lBlaes. Lee. & Socratesll2002l ; 
ISetdl20I2h . Even for m2 ~ mi, the characteristic frequency of 
Kozai process is O(a^ni), while the IPN precession frequency of 
the inner EMRI is ~ 3pni. When a system encounters our reso- 
nance, we have a'^^^ = 3p ^ 1 (see Eg. ill It) and the Kozai pro- 
cess is suppressed by the IPN precession effect (see also Naoz et 
al. 2012). Note that the semi major axes of the systems shown here 
are not constant of motion, and the orbital averaging associate with 
the Kozai mechanism cannot be applied here. These systems lay 
below the stability criterion presented in Lithwick & Naoz 201 1. 

As for the coplanar orbits, we present the capture rates 
(Fig.I2) and the gaps of the eccentricities (Fig.l3). The analytical 
predictions based on the Hamiltonian approach show good agree- 
ments with the numerical results at e < 0.4. 

The favourable results in this section could be regarded as ad- 
ditional supports for validity of our Hamiltonian approach extended 
for the relativistic resonance. 



9 SUMMARY 

We have studied dynamics of the resonant state A2 — ^1 ~ const 
for a triple system composed by an EMRI (CO-l-MBH) and an ad- 
ditional outer MBH. This resonant state is supported by the rela- 
tivistic apsidal precession of the inner EMRI, and does no depend 
on its mean anomaly Ai . As a result, the two orbits can become 
hierarchical with a = axja^ <^ 1, and then the two masses of 
the MBHs mo and m2 can become comparable, in contrast to the 
standard MMRs where we have mi, m2 ^ mo due to dynamical 
stability of orbits with a = 0(1). 
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Figure 12. Resonant capture rate of an inner EMRI with the inclination 
angle i = 60° . The solid curves are given by the Hamiltonian approach 
with the second order mode. The circles are results from numerical results 
(each from 20 runs) with the minimum duration 02/10 of the resonant state, 
and the additional open squares are for 02/20 (shown only if different). 
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Figure 13. Correspondences of the eccentricities around the resonant en- 
counters. The results are given for models I and II, as in Fig. 8, but now for 
inclined orbits with i = 60° . 



As a preliminary analysis, in §4, we discussed the dominant 
order of the resonant interaction for our state A2 — roi ~ const. 
Due to the orbital hierarchy, dependencies on the parameter a{<^ 
1) play a critical role to assess the relevant terms, and the second- 
order one oc Q^e^ (e; the eccentricity of the EMRI) could become 
more important than the first-order one oc ea"^. This result is re- 
markably different from the standard MMRs for which the param- 
eter Q = 0(1) is less important. 

In §5, we derive the mapping from the resonant triple sys- 
tems to the rescaled one-dimensional Hamiltonian for the state 
A2 — TUi ~ const. We basically followed the framework of the 
Hamiltonian approach explained in the literature, but payed special 
attention to the term associated with the relativistic apsidal preces- 
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sion zuir- Here the dependence -cjir oc (1 — e^)~^ on the eccen- 
tricity e is the key element for the structure of the derived Hamilto- 
nian. The mapping from the original EMRI-MBH triple system to 
the rescaled Hamiltonian becomes regular even with the test parti- 
cle limit mi — where the difference from the standard MMRs 
would become clear. 

Then, based on the derived mapping, we made analytical pre- 
dictions on the dynamical evolution of the state A2 — rui ~ const 
around the resonant encounter and compare them with numerical 
simulations. 

In §6, we studied the resonant capture rate as a function of the 
eccentricity e. For the analytical rate, we incorporated the mapping 
derived in §5 with the expressions given by Borderies & Goldreich 
(1984) for the rescaled Hamiltonian. We found that our analytical 
rates show reasonable agreements with numerical results for eccen- 
tricity e < 0.4 where we can perturbatively deal with the effects of 
the eccentricity e. 

In §7, we studied the gap of the inner eccentricity when the 
capture is failed. With the rescaled Hamiltonian, this characteris- 
tic phenomena can be understood as a sudden change of periodic 
motion at a separatrix crossing. We showed that our analytical pre- 
dictions matches numerical results well. 

Finally, in §8, we discussed relatively inclined orbits. By eval- 
uating dependence of the disturbing function on the inclination an- 
gle, we can derive the relevant expressions required for the map- 
ping between the inclined triple system to the rescaled Hamilto- 
nian. Again, our analytical prediction reproduces numerical results 
well for e < 0.4. 

In this paper, setting EMRI-MBH triple systems as our con- 
crete astrophysical targets, we studied the hierarchical resonant 
state A2 — TUi ~ const induced by relativistic apsidal preces- 
sion. Similar analyses might be useful for purely Newtonian sys- 
tems such as a planet orbiting around one component of binary 
stars. Also in these cases, the mapping could be well behaved in 
the test particle limit mi — 0, due to preferred dependencies of 
the apsidal precession rates on the inner eccentricities. 

This work was supported by JSPS (20740151, 24540269) and 
MEXT (24103006). 
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APPENDIX A: FIRST ORDER RESONANCE 

In Eq.JI9l) we expand the resonant potential ( up to second-order 
in eccentricity e. In §5, we derive the mapping to the simplified 
Hamiltonian only including the second order mode cos 277. Here 
we present formulae keeping the first-order mode oc e cos rj alone 
in the test particle limit mi — 0. The basic procedure is the essen- 
tially the same as §5. After some algebra, the rescaled Hamiltonian 
is given in the forms 



with (j> = r]. The momentum $ is related to eccentricity as 

$ = Ee^ 

with the coefficient 

24/3^10/9 



E 



316/952/3^2/3^4/3 ■ 



(Al) 



(A2) 



(A3) 



Meanwhile the transition speed dS/dr is expressed with the origi- 
nal parameters of the triple system as 



d5_ 
dr 



(A4) 



Note that Eqs.( IA3t and ( IA4b are regular in the test particle limit, as 
in the case for the second-order mode. 



APPENDIX B: ANALYTICAL FORMULAE FOR THE 
RESONANT CAPTURE RATES 

Borderies & Goldreich (1984) analytically evaluated the resonant 
capture rate for a simplified Hamiltonian 
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12 A^. Seto 



H' = $' + (5$+(-l)'=-l' 



fe/2 



(Bl) 



with fc = 1 (first-order resonance) and fc = 2 (second-order reso- 
nance). Here a capture means that the motion of the angle shifts 
from a rotation in the full angular range [0, 2tt\ to a libration within 
a limited range. We should note that, as partially demonstrated in 
§5, perturbative expansion for the eccentricity e is crucial to derive 
the rather simple forms (IBIK 

The goal in this appendix is to provide their capture rates as 
functions of the initial momentum We do not intend to fully 
explain the arguments in Borderies & Goldreich (1984), but con- 
cisely presents their main results. We basically deal with the two 
cases k — 2 and k — 1 separately in §B1 and B2 below. Before 
going into analyses specific to each mode, we firstly describe the 
features common to fc = 2 and 1. Here we omit the label fc, if 
unnecessary 

We consider a case when the parameter S decreases adiabati- 
cally from <5 3> 1 down to 5 <^ —1. At an initial epoch 5i„ 2> 1, 
far from resonance, the structure of a curve H' = const is simple. 
We have $ ~ "l?i„ = const with the conjugate variable (j) rotating 
in the full range [0, 2n]. When the parameter 5 becomes less than a 
critical value Scr (depending on the parameter fc), the Hamiltonian 
H' now has an inner separatrix, in addition to an outer one. We 
denote the areas inside these two separatrixes by 2n^i (inner) and 
27r$o (outer), as functions of the epoch S. For the critical value 
5 = Scr, we put t&cr = ^o{Scr) (with the identity $/(5cr) = 0). 

From conservation of an adiabatic invariant 
jLandau & Lifshitzl Il969l) . the system with an initial value 
$in < $cr is captured into the resonance at 100%, namely 
with the capture rate of P($in) = 1. In terms of the decreasing 
parameter S, this happens at 5 > Scr- 

But, for ^ ^cr, the resonant encounter occurs at 5 ^5 Scr 
with the capture rate Pfc($i„) < 1 that was estimated using an ar- 
gument based on energy balance (see e.g. Goldreich & Peale 1966). 
If the capture failed, the system starts to rotate near the inner sepa- 
ratrix and soon relaxes to a simple rotation state with the magnitude 
$ = $/ (<!>/: evaluated at the resonant encounter). 

Below, for fc = 1 and 2, we separately provide the formulae 
that were given in Borderies & Goldreich (1984) but appropriately 
adjusted for the specific forms of our Hamiltonian ( IB It . 



Bl Second-Order Resonances; fc = 2 

The critical values are Scr = — 1 and "I>cr ~ 1. The capture rate for 

^in > ^cr is 



1 + 



^ 2arcsin[{-i)-l/2] 

where the parameter S is related to as 

^in = ^oiS) 

with the explicit form $0 {S) as 

i-S-iy/^ S r-K 



(B2) 



(B3) 



■i>o{S) 



^||+arcsin[(-5)-^/^]}. (B4) 



The expression for $/ (S) is given in a similar complicated form. 
But, using a formula of a definite integral, we can obtain a simple 
relation (Malhotra 1988) 

$0 + = -5. (B5) 
We plot the rate for fc = 2 in Fig IB II with the solid curve. 
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Figure Bl. Resonant capture rates for the scaled Hamiltonian (Bl) with 
fc = 1 (first-order resonance) and k = 2 (second-order resonance). 



B2 First-Order Resonances; fc = 1 

The critical values are Scr ~ —3/2 and $cr ~ 3/2. The rate for 
> $cr is given as 

2 



Pl('l>,n) 



(B6) 



2arcsin[(-2«u/3)--'*/2] 

where the two parameters S and u are the solutions of the following 
two equations 

-3/2] 

O- I 7T I — Ann \ 

27r$„i 



— < — + arcsm 



9_ 

2Su 

27 „ 2 

= 6 — u . 



-2Su 



3 

-| 1/2 



4:S^U 

In Fig lBll the capture rate for fc = 1 is shown with the dashed 
curve. 



APPENDIX C: ADDITIONAL CANONICAL 
TRANSFORMATION 

The rescaled Hamiltonians H in Eqs.l|49ll and l lAlb are given for 
the variable (j> (closely related to 77) and its conjugate momentum 
<I>(oc e^). The corresponding canonical equations are written as 

dtp _dH_ 9£ _ dH 

In the literature, an additional canonical transformation {(p, $) 
(a;, y) is frequently introduced as follows 



(CI) 



X = V2$ cos (f), y — \/2$ sin (; 



(C2) 



In the new phase space, the angle between the x-axis and the point 
(x, y) is identical to the original variable (j), and the distance to the 
origin is proportional to the eccentricity as y'a;^ 4- y2 q<- 

oc 

e. With these new variables, the rescaled Hamiltonians H become 
simple polynomials as 

'2|2\2 2,2 2 2 

H=\^-±y-\ +s^_4^ + ^_]L (C3) 



(C4) 



for the second-order one ( 149) , and 

' 2 , 2 \ 2 2 I 2 



V2 



© 0000 RAS, MNRAS 000, 000-000 
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for the first-order one (lAll . 

The canonical variables [x, y) have geometrically intuitive 
meanings, and, in many cases, they are more convenient to an- 
alyze the dynamics of the rescaled Hamiltonians themselves, as 
widely done in preceding studies (Borderies & Goldreich 1984; 
Peale 1986; MD). But we concentrate on the issues more related to 
the mapping between the Hamiltonians and the EMRI-MBH triple 
systems, and stay away from the new variables [x, y) in the rest of 
this paper. 
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